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Abstract 

We study dynamics of the equation of state parameter for a dark energy compo- 
nent non-minimally coupled to induced gravity on a warped DGP brane. We show 
that there are appropriate domains of the model parameters space that account 
for crossing of the phantom divide line. This crossing, which is possible for both 
branches of the scenario, depends explicitly on the values of the non-minimal cou- 
pling and warp factor. The effect of warp factor appears in the value of the redshift 
parameter at which phantom divide line crossing occurs. 
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1 Introduction 



Recent observational data from CMB temperature fluctuations spectrum, Supernova type 
la redshift-distance surveys and other data sources, have shown that the universe is cur- 
rently in a positively accelerated phase of expansion and its spatial geometry is nearly fiat 
[1]. Nevertheless, there is not enough standard matter density in the universe to support 
this flatness and accelerated expansion. Therefore, we need either additional cosmological 
components or modify general relativity at cosmological scales to explain these achieve- 
ments [2]. Multi-component dark energy with at least one non-canonical phantom field^ is 
a possible candidate of the first alternative. This viewpoint has been studied extensively 
in literature ( see [5,6] and references therein ). There are some datascts (such as the 
Gold dataset) that show a mild trend for crossing of the phantom divide line by equation 
of state (EoS) parameter of dark component. The equation of state parameter in these 
scenarios crosses the phantom divide line ( = ^ = — 1) at recent redshifts and current 
accelerated expansion requires a; < — |. In fact, recent observational data restrict u!{z) to 
be larger than —1 in the past and less than —1 today. The current best fit value of the 
equation of state parameter, using WMAP five year data combined with measurements of 
Type la supernovae and Baryon Acoustic Oscillations in the galaxy distribution, is given 
by —0.11 < 1 -I- a; < 0.14 ( with 95 percent CL uncertainties) [7]. It is accepted that 
crossing of the phantom divide fine occurs at recent epoch with z ~ 0.25 [5,6], although 
this value is model dependent. Currently, models of phantom divide line crossing are so 
important that they can realize that which model is better than the others to describe the 
nature of dark energy [6] . Although this crossing cannot be explained just by one scalar 
field [8], generalization to multi-field case or non- minimal coupling with gravity provide 
enough space to achieve such a crossing [3,5,6]. Lorentz invariance violating fields are 
other alternative dark energy components with capability to cross the phantom divide 
line by EoS parameter in a fascinating manner [9]. Also it is possible to reconstruct a 
scalar-tensor theory of gravity in an accelerating universe where a phantom behavior can 
realized [10]. The cosmological constant ( with u!{z) — —1) is the simplest candidate 

^It is important to note that phantom fields are not consistent due to violation of the null energy con- 
dition and instabilities ( see for instance [3]). However, theoretically they provide a good candidate with 
negative pressure to realize late-time accelerated expansion. Recently, it has been shown that phantom- 
like behavior can be realized without introducing any phantom matter in some specific braneworld models 
[4]. We note also that in our model due to its wider parameter space, it is expected essentially that the 
null energy condition is fulfill in at least some subspaces of the model parameter space. 
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for dark energy [11,12]. However this scenario suffers from some difficulties sucli as lack 
of physical motivation, huge amount of fine-tuning to explain cosmological accelerated 
expansion and no dynamics for its equation of state [12]. So it seems worthwhile to probe 
alternative dynamical models. 

Another alternative to explain current accelerated expansion of the universe is exten- 
sion of the general relativity to more general theories on cosmological scales ( see [13] and 
references therein). DGP ( Dvali-Gabadadze-Porrati) braneworld scenario as an infra-red 
(IR) modification of general relativity, explains accelerated expansion of the universe in 
its self-accelerating branch via leakage of gravity to extra dimension [14]. In this model, 
the late-time acceleration of the universe is driven by the manifestation of the excruciat- 
ingly slow leakage of gravity off our four-dimensional world into an extra dimension [15]. 
In this scenario the EoS parameter of dark energy never crosses the uj{z) = —1 line, and 
universe eventually turns out to be de Sitter phase. However, in this setup if we use a 
single scalar field (ordinary or phantom) on the brane, we can show that EoS parameter 
of dark energy can cross the phantom divide line [16]. It has been shown that DGP model 
with a quintom dark energy fluid ( a combination of quintessence and phantom fields in 
a joint model) in the bulk or brane, accounts for accelerated expansion and the phantom 
divide line crossing [17]. In a braneworld setup with induced gravity embedded in a bulk 
with arbitrary matter content, the transition from a period of domination of the matter 
energy density by nonrelativistic brane matter to domination by the generalized dark ra- 
diation provides a crossing of the phantom divide line [18]. In this setup there is no need 
to introduce additional scalar field as a dark component. Recently, phantom-like behavior 
in a brane- world setup with induced gravity and also curvature effects have been reported 
[19]. On the other hand, Gauss-Bonnent braneworld scenario with induced gravity dose 
not need introducing any scalar field to account for this crossing. In other words, the 
combination of the effect of Gauss-Bonnent term in the bulk and induced gravity term 
on the brane behaves as dark energy on the brane [20]. Quintessential scheme can also 
be achieved in a geometrical way in higher order theories of gravity [21]. 

Crossing of the phantom divide line by a minimally coupled scalar field on the DGP 
braneworld has been studied by Zhang and Zhu [16]. In this setup, there are two possible 
cases: for ordinary scalar field EoS of dark energy crosses from a;>— ltoa;<— lin 
normal ( non self-accelerating) branch and for phantom field EoS of dark energy crosses 
from a;<— ltoa;>— lin self-accelerating branch of DGP scenario. As an important 
generahzation, the Randall-Sundram II model [22] combined with DGP scenario provides 
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a rich structure which has been called warped DGP braneworld in literature [23] . In this 
model, an induced curvature term appears on the brane in the RS II model. 

Our motivation in this paper is to consider a scalar field non-minimally coupled to 
induced gravity on the warped DGP braneworld as a dark energy component and inves- 
tigate the roles played by non-minimal coupling and the warp effect in the dynamics of 
the equation of state parameter. We study dynamics of the equation of state parameter 
focusing on the crossing of the phantom divide line in this setup. We show that this 
crossing is possible for a suitable range of the model parameters and especially for some 
specific values of the non-minimal coupling and warp factor. More specifically, we show 
that for ordinary scalar (quintessence) field in the self-accelerating branch of the warped 
DGP model and with positive values of the non-minimal coupling, the EoS of dark energy 
runs from below —1 {ujde < —1) to above —1. In normal branch of this warped DGP 
model with positive and negative values of the non-minimal coupling, the EoS parameter 
runs from above —1 to below —1 which is supported by recent observations. For phantom 
field we have crossing of the phantom divide line in both branches of the warped DGP 
setup and in both of these cases the EoS parameter runs from above —1 to below — 1. 
But, in self- accelerating branch of the model we have crossing behavior with positive and 
negative values of the non- minimal coupling in normal branch and crossing occurs just 
with negative values of the non-minimal coupling. 

2 A Dark Energy Model on the Warped DGP Brane 
2.1 Warped DGP Brane 

Let us start with the action of the warped DGP model as follows 

(2) 



iZKk J J brane 



bulk -2/t5 



Here S^uik is the action of the bulk, S^r-ane is the action of the brane and S is the total 
action. with A = 0, 1, 2, 3, 5 arc coordinates in the bulk, while x^^ with = 0, 1, 2, 3 
are induced coordinates on the brane. is the 5-dimensional gravitational constant. ^^^R 
and ^^^Cm are 5-dimensional Ricci scalar and matter Lagrangian respectively. is trace 
of the extrinsic curvature on cither side of the brane. ChraneiQaiSii^) is the effective 4- 
dimensional Lagrangian. The action S is actually a combination of the Randall- Sundrum 
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II and DGP model. In other words, an induced curvature term is appeared on the brane 
in the Randall- Sundrum II model, hence the name warped DGP Braneworld [23] . Now 
consider the brane Lagrangian as follows 



(3) 



where // is a mass parameter, R is Ricci scalar of the brane, A is tension of the brane 
and Lm is Lagrangian of the other matters localized on the brane. Assume that bulk 
contains only a cosmological constant, ^^^A. With these choices, action (1) gives either a 
generahzed DGP or a generalized RS II model: it gives DGP model if A = and ^^^A = 0, 
and gives RS II model ii fj, — [23] . The generahzed Priedmann equation on the brane is 
as follows [23] 



rr2 k 1 
^ + - = ^ 



p + Pofl + a 



(4) 



where e = ±1 is corresponding to two possible branches of the solutions in this warped 



DGP model and A 



1/2 



where An 



1 



1/2 



V 



6m§ 



POM" 

manifold and 



with < ?7 < 1 which is related to the warped geometry of the bulk 
Pq = m\ + 6^. By definition, mx — X^^'^ and — k^^^^. Sq is an integration constant 
and corresponding term in the generalized Priedmann equation is called dark radiation 
term. We neglect dark radiation term in which follows. In this case, generalized Priedmann 
equation (4) attains the following form. 



k 



1 



3/x^ 



p + po + epo(^Al + 



2r]p\i/2 



Po 



(5) 



where p in our case is the total energy density, including scalar field and dust matter on 
the brane 

P<p + Pdm- (6) 

In which follows we construct a dark energy model on the warped DGP setup. 



2.2 Quintessence Field 

Now we consider a quintessence scalar field non-minimally coupled to induced gravity on 
the warped DGP brane as a candidate for dark energy. The action of this non-minimally 
coupled scalar field is given by 



brane 



-9 



1 



(7) 
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where ^ is a non-minimal coupling and R is Ricci scalar of the brane. We have assumed 
a conformal couphng of the scalar field and induced gravity. The scalar field will play the 
role of dark-energy component on the brane. Variation of the action with respect to </? 
gives the equation of motion of the scalar field 

dV 

(p + 3i/0 + + — = 0. (8) 
dip 

The energy density and pressure of the non-minimally coupled scalar field are given by 

= + V{^) + ^H^^ + 3ei/V' (9) 

= - y{^) - 2^(W + 2^^<^ + - i^\2H + m^) (10) 

In which follows, by comparing the modified Priedmann equation in the warped DGP 
braneworld with the standard Priedmann equation, we deduce a definition for equation 
of state of dark energy component. This is reasonable since all observed features of dark 
energy are essentially derivable in general relativity [16,20]. The standard Priedmann 
equation in four dimensions is written as 

+ ^ = ^{Pdm + Pde), (11) 

where pdm is the dust matter density, while pde is dark energy density. Comparing this 
equation with equation (5) for a spatially flat universe ( k — 0), we find 

Pde = + Po + epo (aI + 2r]-^) ~\ (12) 

Po^ 

Conservation of the scalar field effective energy density leads to 



dp^ 



+ 3H{p^ + p^)^0. (13) 



dt 

We note that the non-minimal coupling of the scalar field to the Ricci curvature on the 
brane preserves conservation of the scalar field energy density^. 



^Note that authors of Rcf. [24] have treated this conservation in relatively different way. They have 
defined a total energy-momentum tensor consist of two parts: a pure (canonical) scalar field energy- 
momentum tensor and a non-minimal coupling-dependent part. The total energy density defined in this 
manner is then conserved. In our case, we have included all possible terms in equations (9) and (10) 
from beginning and evidently total energy density defined in this manner is conserved too. In fact, it 
is simple to show that our equations (9) and (10) arc equivalent to p*"* and P*°* of Ref. [24] if we set 
«(</?) = i(l — ^(p^) ( see also [25] for more detailed discussion). 
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Since the dust matter obeys the continuity equation and the Bianchi identity keeps 
vahd, dark energy itself satisfies the continuity equation 



dp, 



'de 



dt 



(14) 



where pde denotes the pressure of the dark energy. Note that this is just a definition and 
by assuming vahdity of this definition we can obtain effective pressure of dark energy as 
well as an effective equation of state. Now, the equation of state for this dark energy 
component can be written as follows 



Pde 1 1 dlnpde 
Pde 3dln(l + ^) 



(15) 



where by using (13) and (14) we find 



dlnpde 
d\n{l + z) 



p^+p^ + erj(Al + 2r]P^^P'"' 



Po 



P^+P^ + Pdm 



(16) 



There are three possible alternatives in this setup: if | Ji!^(i+z) > 0' have a quintessence 
model; if | din{i+z) ^ ^^'^ model is phantom and if | rftnfl+z) ~ ^' dark component is 
a cosmological constant. Evidently, in this setup non-minimal coupling of the scalar field 
and induced gravity plays a crucial role supporting or preventing phantom divide line 
crossing. In this respect, the differences between the minimal and non- minimal setups 
will be more clear if we write the explicit dynamics of the equation of state parameter. We 
also discuss the effect of warp factor on the dynamics of the equation of state parameter 
in forthcoming arguments. We choose the following exponential potential with motivation 
that this type of potential can be solved exactly in the standard model 



y = yoexp(-A^), 
A* 



(17) 



where Vq, A and p, are constant. 

Differentiation of the logarithm of dark energy effective density with respect to ln{l + z) 
yields 



din Pde _ _3_ 
din (1 + 2;) Pde 



er](^Al + 2ri^ 



if' + V{ip) + G^Hifip + S^H'ip' + Pdrn\-\ 



Po 



(18) 
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To study the behavior of the EoS parameter of dark energy component, we need to the 
exphcit form of (p in terms of other quantities which can be deduced from equation of 
motion as given by (8). On the other hand, Priedmann equation given by (5) now takes 
the following form 



{ii^ + gfH^ + 2/(3/^2 + g)H'' + f- 2/(3/^2 + g) + 2rjpog 



where 



+ ( - 2// + povf)H - 2r]po{l - po) - Po^o + = (19) 

I = ^<yf^ + V{ip) + pdm + Po, 



and 



Equation (20) is a quadratic equation in terms of and in principle has four roots for 
H . We show these roots as /ii, /i2, /is and /14. After numerical calculation, we found that 
two of these roots, say, hi and /i2 are unphysical and excluded by observational data. The 
other two roots, /is and /),4, are physical solutions corresponding to the generalized normal 
branch (with e = — 1) and the self-accelerating one (with e = — 1). As we will show these 
solutions have the capability to account for phantom divide line crossing. We introduce 
a new parameter defined as s = — ln(l + z) and rewrite dark energy equation of state 
parameter as follows 

Idlnpde , 
"^^^ = -'-3^1^- 

Now, we analyze the behavior of Wde versus s to investigate cosmological implications 
of this scenario. Using equation (18), we see that in the minimal case ( with ^ = ) 
and neglecting warp effect, if we choose the sign of e to be negative, remaining terms 
have suitable combination of signs so that it is possible to cross the phantom divide line 
by the EoS parameter. In the non-minimal case, however, it is not simple to conclude 
that there is crossing of the phantom divide line or not just by defining the sign of e, 
since in this case non-minimal coupling itself has a crucial role in the dynamics of the 
equation of state parameter. We consider ^ as a fine-tuning parameter in this case. It 
is important to note that in the absence of the scalar field on the DGP setup, there is 
no crossing of the phantom divide line on the self-accelerating branch even in the warped 
DGP scenario. Nevertheless, normal branch accounts for crossing of the phantom divide 
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line in this situation. In which follows, to calculate EoS parameter u!{z), we define some 
dimensionless density parameters such as ( see the paper by Sahni and Shtanov in Ref. 
[4] and also Ref. [15]) 

where we have assumed that p is given by the sum of the energy densities of different 
components labeled by a with a constant EoS parameter, Ua- Also we define 

A; 1 A (^^A 

H^ai 3p,^H^ 6H^ 

where Tc is DGP crossover distance. In this case, Priedmann equation can be written as 
follows 



= ^^(1 + Zf + ildmiil + Zf - 1] + ^4(1 + 0)3(1+") - 1] + l^A + 



Hi 



±2 



VLr^ [VLdm[{l + Zf - 1] + + 0)3(1+") - 1] + + + 



1/2 



(23) 



where ± stands for two possible embedding of the brane in the bulk. The constraint 
equations for cosmological parameters are given by 

2 



1 — VLu -\- Vlf 



(24) 



We can define also Q,ki as the present value of the scalar field kinetic energy density |</?^ 
over the critical density defined as pc = Si^Q ( with SttG = 1). Note that Vt^i and the 
non- minimal coupling parameter are hidden in the definition of the Vl^. If we change 
the values of these parameters in appropriate manner, the redshift at which crossing of 
the phantom divide line occurs will change since it is a model dependent quantity in this 
respect. In table 1, we have obtained some reliable ranges of the non-minimal coupling 
to have crossing of the phantom divide line in this setup. Observational data show that 
crossing of the phantom divide line is occurred in redshift ^ ~ 0.25, so wc have obtained 
the values of ^ which are correspond to this redshift in the last column of the table 1. We 
have not excluded the negative values of the non-minimal coupling from our analysis. In 
fact, these negative values arc theoretically interesting, corresponding to anti-gravitation 
( see [25] for further discussion). The results of the numerical calculations are shown in 



''Note that this is a model dependent value but this value is suitable for our purposes in forthcoming 
arguments. 
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Table 1: Acceptable range of ^ to have crossing of the phantom divide line with 
quintessence field ( constraint by the age of the universe). 



e 




Acceptable range of ^ 


The value of ^ for z=0.25 


+ 1 


negative 


-0.605 < C < 


-0.438 


+ 1 


positive 


no crossing 




-1 


negative 


-0.83 < C < 


-0.522 


-1 


positive 


< ^ < 0.148 


0.124 



figures 1, 2, 3 and 4 for two branches of this DGP-inspired model and with different values 
of the non-minimal coupling ^. In this figures, the best ranges of the values for ^ to have 
a reliable model in comparison with observational data are obtained. Note that in all of 
our numerical calculations we have assumed fij^j = 0.01, fi^c = 0.01, Vt^ = 0.3, Vt^ = 0, 

= = 0.7, Qjy = l, Ao = l, Ho = l, 11 = 1 and r] = 0.99. 
Figure 1(a) shows that for self-accelerating branch of the model ( with e = +1) and with 
negative values of the non-minimal coupling, there is a crossing of the phantom divide 
line by the equation of state parameter( note that self-accelerating branch of this DGP- 
inspired model accounts for late-time accelerated expansion). As we know, the EoS of a 
single scalar field in standard relativity never crosses the phantom divide line ( see for 
instance the paper by Vikman in Ref. [8]). Also, in DGP scenario with a canonical scalar 
field on the brane, there is no crossing of the phantom divide line in the self- accelerating 
branch of the model ( sec [16] and [27]), but in our model we have such a crossing due to 
the existence of the non- minimal coupling and warp effect. In our case, if we choose for 
instance ^ = —0.438, we have a phantom divide line crossing at the point with s = —0.22 
corresponding to z = 0.25 in agreement with observations. However, this crossing occurs 
for negative values of the non-minimal coupling. On the other hand, here the EoS runs 
from below —1 to above — 1 ( from phantom to quintessence phase) and therefore avoids 
big-rip singularity. 

The effect of the warp factor on the dynamics of the dark energy component can be 
explained by the variation of rj parameter as shown in figure 1(6). As this figure shows, for 
sufficiently small values of rj, equation of state parameter, uj, crosses the phantom divide 
line in relatively small values of redshift. In figure 1(c) we plotted cude 

for DGP^+^ branch 
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Figure 1: a) In the self-accelerating branch of the model and with negative values of the 
nonminimal coupling, the EoS parameter crosses the phantom divide line. For instance, with 
^ = —0.438 this crossing occurs at s = —0.22 or z = 0.25. b)The role played by r] ( which is 
related to warp effect) on the crossing of the phantom divide line. For sufficiently small values 
of r], equation of state parameter, oj, crosses the phantom divide line in relatively small values 
of redshift. For example, with ^ = —0.438, the EoS of dark energy crosses lo = —1 line with 
T] = 0.99 at s ~ —0.22 or z = 0.25, while for rj = 0.50 this crossing occurs at s ~ —0.212 or 
z = 0.236 and for rj = 0.10 this crossing occurs at s ~ —0.198 or z = 0.218. c) Equation of 
state parameter, w, versus s and tj with = —0.438 in a three dimensional plot and within self- 
accelerating branch. In self-accelerating branch with negative non-minimal coupling, crossing of 
the phantom divide line occurs from phantom to quintessence phase. 
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of the model with ^ = —0.438 with respect to the parameters s and r). In this figure, r) is 
restricted to the interval < 77 < 1. 

In figure 2 we see that for self-accelerating branch of the model and with positive 
values of the non-minimal coupling, there is no crossing of the phantom divide line. 

Figure 3(a) shows that for normal branch ( with e = — 1) of the model and with 
negative values of ^, there is crossing of the phantom divide fine also. For instance, with 
^ — —0.522, we have a crossing at the point with s — —0.22 corresponding to ^ = 0.25 in 
agreement with observational data. The other note is that EoS transits from Ude > — 1 to 
iOde < — 1 which is supported by recent observation. Figure 3(6) shows that by decreasing 
the values of r], equation of state parameter, 00, crosses the phantom divide line in relatively 
small values of redshift. In figure 3(c) we plotted cude for DGP^~^ branch with ^ — —0.522 
and with respect to the parameters s and rj. 

Figure 4(a) shows that EoS parameter of dark energy crosses the phantom divide line 
in the normal branch (e = —1) of the model with positive ^. For ^ = 0.124 we have a 
crossing at the point with s = —0.22 corresponding to z = 0.25. Here, the EoS of dark 
energy transits from Ude > — 1 to ujde < — 1 in agreement with the recent observation 
which show crossing from quintessence to phantom phase. Figure 4{b) shows that by 
reduction of the values of i], uj crosses the phantom divide line in relatively lower values 
of redshift. In figure 4(c) again we plotted Ude for DGP^~^ branch with ^ = 0.124 with 
respect to the parameters s and 77. 
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Figure 2: With a quintessence field on the warped DGP brane, there is no crossing of the 
phantom divide fine in self-accelerating branch of the model with positive values of the non- 
minimal coupling. This means that inclusion of the warp factor cannot produce crossing of the 
phantom divide in this case. 
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C) E = -1 , ^ =-0.522 




Figure 3: a) In the normal branch of the model, there is a crossing of the phantom divide line 
with negative values of the nonminimal coupling. This crossing occurs from quintessence to 
phantom phase. For instance, the EoS parameter crosses the oj = —1 line for ^ = —0.522 at 
s = —0.22 or z = 0.25. b)The role played by r/ factor. For sufficiently small values of r/, equation 
of state parameter, a;, crosses the phantom divide line in relatively small values of redshift. 
For example, with ^ = —0.522, the EoS of dark energy crosses lo = —1 line with r] = 0.99 at 
s ~ —0.22 or z = 0.25, while for r] = 0.7 this crossing occurs at s ~ —0.203 or z = 0.225 and 
for T] = 0.5 occurs at s —0.189 or z = 0.208. c)Equation of state parameter w, versus s and rj 
with ^ = —0.522 in a three-dimensional plot and for normal branch of the model. Note that in 
this case, crossing runs from quintessence to phantom phase. 
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c)e = -1 ,5=0.124 




Figure 4: The same as figure 3 but now with positive values of the non-minimal coupling. a)In 
this case crossing of the phantom divide line runs from quintessence to phantom phase. For 
^ = 0.124 this crossing occurs at s = —0.22 or z = 0.25. b)The role played by rj on the crossing 
of the phantom divide line. For sufficiently small values of r], equation of state parameter crosses 
the phantom divide line in relatively small values of redshift. For instance, if ^ = 0.124, the EoS 
of dark energy crosses to = —1 line for rj = 0.99 at s ~ —0.22, or z ~ 0.25, while for rj = 0.7 this 
crossing occurs at s « —0.16 or z = 0.173 and for rj = 0.5 occurs at s ~ —0.10 or z ~ 0.105. 
c)Equation of state parameter, w, versus s and rj with ^ = 0.124 in a 3-dimensional plot. 
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2.3 Phantom field 



Now we investigate dynamics of a phantom field non-minimally coupled to induced gravity 
on the warped DGP brane. Most of the techniques and discussions for this case are similar 
to the previous subsection. The action of the model is 



brane 



(25) 



where a is a phantom field. Variation of the action with respect to a gives the equation 
of motion of the phantom field 



dV 

a + 3H& - $Ra - — = 0. 

da 

The energy density and pressure of this phantom field are given as 

1 



and 



Pa = -^CT^ + V{a) + 6^Haa + S^H^a^ 



Pa = --(7^ - V{a) - 2^{aa + 2aHa + a') - ^a'{2H + 3H'). 



(26) 



(27) 



(28) 



To compare with the results corresponding to the quintessence field, we assume the same 
type of potential 

y((7) = yoexp(-A-) (29) 

n 

where Vq, A and /i are constant. Differentiation of the effective energy density of phantom 
field with respect to ln{l + z) is given by 



dp, 



de 



d]n(l + z) 



3 

Pde 



^ Po 



(30) 



where a can be deduced from equation of motion of a, (26). On the other hand, Priedmann 
equation now takes the following form 



ip' + gYH^ + 2f{3p' + g)H' + f- 2l{3p' + g) + 2rjpog 



+ {-2fl + povf)H - 2r,po{l - po) - plAl + = 



(31) 



where by definition 



g = -3iH'a\ 
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Table 2: Acceptable range of $, to have crossing of the phantom divide hne with just one 
phantom field ( constraint by the age of the universe) . 



e 




Acceptable range of ^ 


The value of ^ for z=0.25 


+ 1 


negative 


-0.485 < C < 


-0.366 


+ 1 


positive 


0.055 < ,^ < 0.170 


0.088 


-1 


negative 


no crossing 




-1 


positive 


< ^ < 0.166 


0.166 



I = "2^^ + + Pdm + PO, 

and 

/ = -6^Ha&. 

Similar to the last subsection, there is a forth order equation for H and in principle this 
equation has four roots. Two of these roots are un-physical but two remaining solutions 
are physical and corresponding to two branches of solutions in this DGP-inspired model. 
In table 2, we have obtained some acceptable ranges of non- minimal coupling to have 
crossing of the phantom divide line in this setup and constraint by the age of the universe 
( that is, we have assumed that the age of the universe is ISGyr). 

Figure 5(a) shows that for self-accelerating branch of the model ( with e — +1) and 
with negative values of the non-minimal coupling, there is a crossing of the phantom di- 
vide line by the equation of state parameter. For ^ = —0.366, we have a crossing at the 
point with s = —0.22 corresponding to z = 0.25 in agreement with observations. Here 
the EoS runs from above —1 to below —1 ( from quintessence to phantom phase). Figure 
5(6) shows that for sufficiently small values of r), equation of state parameter, uj, crosses 
the phantom divide line in relatively large values of redshift. In figure 5(c) we plotted oUde 
for DGP^^^ branch with ^ = —0.366 and with respect to the parameters s and rj. 

Figure 6(a) shows that for self- accelerating branch of this model and with positive 
values of ^, there is a crossing of the phantom divide line too. For example, with ^ = 0.088, 
we have a crossing at the point with s = —0.22 corresponding to z = 0.25 in agreement 
with observational data. As another important point, the EoS parameter transits from 
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a) £ = +1 , ^ < b) E = +1 , 5 =-0.366 




0.5 



s 1 



Figure 5: Dynamics of the equation of state parameter with a phantom field on the self- 
accelerating branch of the model, a) With negative values of the nonminimal coupling, the EoS 
of dark energy crosses u = —1 line running from quintessence to phantom phase. For instance, 
with = —0.366, this crossing occurs at s = —0.22 or z = 0.25. b)The role played by rj in 
equation of state of phantom field on the self-accelerating branch of the model. For sufficiently 
small values of r], equation of state parameter crosses the phantom divide line in relatively 
large values of redshift ( in contrast to the case with quintessence field). For example, with 
= —0.366, the EoS of dark energy crosses uj = —1 line with rj = 0.99 at s ~ —0.22 oi z = 0.25, 
while for r] = 0.50 this crossing occurs at s ~ —0.250 or z ~ 0.280 and for i] = 0.10 this occurs 
at s « —0.297 or z = 0.345. c) Equation of state parameter, lo, versus s and r] with ^ = —0.366 
in self-accelerating branch of the model in a 3-dimensional plot. In the self- accelerating branch 
with phantom field, uj = —1 line crossing runs from quintessence to phantom phase. 
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cude > — 1 to cude < — 1- Figure 6(6) shows that by decreasing the values of rj, equation of 
state parameter, cu, crosses the phantom divide hne in relatively large values of redshift. 
In figure 6(c) we plotted LOde for DGP^^^ branch with ^ = 0.088 and with respect to the 
parameters s and 77. In figure 7 we see that for normal or non self- accelerating branch of 
the model and with negative values of the non-minimal coupling, there is no crossing of 
the phantom divide line. 
Figure 8(a) shows the phantom divide line crossing of EoS parameter for normal branch 
(e — —1) of the model with positive values of ^. Here, the EoS of dark energy transits 
from ujde > —1 to Ude < —1. Figure 8(6) shows that by reduction of the values of i], u> 
crosses the phantom divide line in relatively smaller values of redshift. In figure 8(c) again, 
we plotted ujde for normal branch of the model with ^ = 0.166 versus the parameters s 
and T). 
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Figure 6: The self-accelerating branch of the model with positive values of the non- minimal 
coupling in the presence of just one phantom field on the brane. The situation is similar to 
previous figure and phantom divide line crossing runs from quintessence to phantom phase, a) 
With positive values of the nonminimal coupling, the EoS of dark energy crosses oj = —1 line for 
^ = 0.088 at s = —0.22 or z = 0.25. b)The role of rj on the crossing of the phantom divide line. 
As previous figure and contrary to quintessence case, for sufficiently small values of r] equation 
of state parameter crosses the phantom divide line in relatively large values of redshift. For 
example, with ^ = 0.088, the EoS of dark energy crosses to = —1 line with r] = 0.99 at s ~ —0.22 
OT z = 0.25, while for rj = 0.5 this crossing occurs at s ~ —0.29 or z = 0.33 and for r/ = 0.1 this 
occurs at s ~ —0.44 or z = 0.55. c) A 3-dimensional plot of the equation of state parameter, uj, 
versus s and r] with ^ = 0.088. 
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Figure 7: With a phantom field on the warped DGP brane, there is no crossing of the phantom 
divide fine in normal ( non self-accelerating ) branch of the model with negative values of the 
non-minimal coupling. Comparison with figure 2 for a quintessence field on the warped DGP 
brane, shows the difi^erences between two situation. 
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c)E = -1 ,^=0.166 




Figure 8: In the normal branch of the model with a phantom field on the brane, crossing of 
the phantom divide line runs from quintessence to phantom phase. This is a general behavior 
for equation of state of phantom field on the brane independent of the signs of non-minimal 
coupling and e. This is supported by observations too. However, as we have shown in figures 1, 
3 and 4, the situation for quintessence field on the warped DGP brane is different and crossing 
of the phantom divide line depends on the signs of the non-minimal coupling and e. That is, 
while with phantom field on the warped DGP brane, crossing of the phantom divide line runs 
always from quintessence to phantom phase, in the case of quintessence field on the brane this 
running of crossing depends on the sign of the non-minimal coupling and e. 
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3 Summary and Conclusion 



An alternative approach to explain current positively accelerated phase of the universe 
expansion is to use a multi-component dark energy with at least one non-canonical phan- 
tom field. The analysis of the properties of dark energy from recent observations mildly 
favor models where ^ = ^ crosses the phantom divide line, u = —1, in the near past. In 
this paper, we have considered a scalar field non-minimally coupled to induced gravity 
on the warped DGP braneworld as a dark energy component and we have investigated 
the roles played by the non-minimal coupling and the warp effect on the dynamics of the 
equation of state parameter. In this respect, wc have studied the dynamics of equation 
of state parameter focusing on the crossing of the phantom divide line in this setup. As 
it is well-known, in the absence of scalar field there is no crossing of the phantom divide 
fine in self- accelerating branch of the DGP setup and this is the case even in the warped 
DGP scenario. However, in the presence of a scalar field ( minimally or non-minimally 
coupled to induced gravity), it is possible to realize this crossing. We have shown that 
this crossing is possible for a suitable range of the model parameters and especially for 
some specific values of the non-minimal coupling and parameter 77 related to the warp 
effect in this DGP-inspired scenario. In the first stage, we have considered a canonical 
(quintessence) scalar field non-minimally coupled to the induced gravity on the warped 
DGP brane. In this case, we have shown that crossing of the cosmological constant fine 
by the EoS parameter of the quintessence field occurs in both self- accelerating and nor- 
mal branches of the model. For self- accelerating branch of the model( with e — +1 ), 
crossing of the cosmological constant line occurs with negative values of the non-minimal 
coupling parameter and this crossing runs from phantom to quintessence phase. There 
is no crossing behavior in the self- accelerating branch of the model with positive values 
of the non-minimal coupling. On the other hand, for normal branch of the model ( with 
e — — 1 ) , the equation of state parameter of dark energy crosses the phantom divide line 
with negative values of the non-minimal coupling as well as its positive values, but the 
crossing behavior is completely different from the former one ( the self-accelerating branch 
). Indeed, in this case the EoS parameter of dark energy crosses the phantom divide line 
in a different direction; from quintessence to phantom phase and this is supported by 
recent observations. By investigating the role played by the parameter r] ( which is re- 
lated to the warp effect) in both branches of the model, we found that decreasing of the 
effect of Tj factor leads to the result that the EoS parameter of dark energy crosses the 
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cosmological constant line in relativity smaller values of redshift. We should stress here 
that negative values of the non-minimal coupling are interesting at least theoretically 
since they show anti-gravitation. However, recent observational constraints on the values 
of the non- minimal coupling favor positivity of this factor ( see for instance [28] and [29]). 

In the next stage, we have considered a phantom field non-minimally coupled to the 
induced gravity on the warped DGP brane. We have shown that in the self- accelerating 
branch of the model, crossing of the phantom divide line by the EoS parameter of dark 
energy occurs with both signs of the non-minimal coupling and this crossing runs always 
from quintessence to phantom phase. Finally, we have shown that with a phantom field on 
the warped DGP brane, there is no crossing of the phantom divide line in the normal ( non 
self-accelerating ) branch of the model with negative values of the non-minimal coupling. 
In the normal branch with phantom field, by considering positive values of the non- 
minimal coupling parameter, the EoS parameter of dark energy crosses the phantom divide 
line from quintessence to phantom phase supported by observations. With a phantom 
field on the self-accelerating branch of the model and for both signs of the non-minimal 
coupling, reduction of the values of 77 leads to phantom divide line crossing in relatively 
larger values of redshift, but in normal branch reduction of r] leads to crossing in smaller 
values of redshift. 

In summary, with positive values of the non-minimal coupling which is physically more 
relevant, we have shown that: in the self-accelerating branch of this warped-DGP setup 
with just one quintessence field, crossing the phantom divide line cannot be realized. 
However, it is possible to realize phantom divide line crossing in the normal branch of 
the model and this crossing runs from phantom to quintessence phase. With just one 
phantom field on the brane, it is possible to realize phantom divide line crossing with 
positive values of the non-minimal coupling in both branches of this DGP-inspired model 
and this crossing runs from quintessence to phantom phase. Although with a phantom 
field on the warped DGP brane, crossing of the phantom divide line runs always from 
quintessence to phantom phase, in the case of quintessence field on the brane this running 
depends on the sign of the non-minimal coupling and e. Finally, we should stress that 
self- accelerating branch of the DGP scenario suffers from ghost instabilities ( see [30,31]). 
Incorporation of new degrees of freedom such as the non-minimal coupling of the scalar 
field and induced gravity and also warp geometry of the bulk provides a wider parameters 
space in our setup and this wider parameter space may provide a suitable basis to treat 
ghost instabihties. 
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